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Abstract 

In this article we present a novel method of hit time and hit position 
reconstruction in long scintillator detectors. We take advantage of the 
fact that for this kind of detectors amplitude and shape of registered 
signals depends strongly on the position where particle hit the detector. 

The reconstruction is based on determination of the degree of similarity 
between measured and averaged signals stored in a library for a set of well- 
defined positions along the scintillator. Preliminary results of validation 
of the introduced method with experimental data obtained by means of 
the double strip prototype of the J-PET detector are presented. 

PACS:87.57.uk, 07.05.Kf 

1 Introduction 

Recently a new concept of large acceptance Jagiellonian PET (J-PET) system 

was proposed [l]-[^. Unlike all the commercial PET devices using inorganic 
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scintillators as radiation detectors [5]-[^ (usually these are the BGO, LSO or 
LYSO crystals) J-PET is based on the polymer scintillators. This technique of¬ 
fers improvement of the Time of Flight (TOE) resolutior0 and also constitutes 
a promising solution for a whole-body PET imaging. In the case of J-PET 
annihilation gamma quanta are registered by means of long scintillator strips 
read out from both ends by photomultipliers. This allows for the determination 
of position and time of the gamma quanta reaction based predominantly on 
the time measurement. Therefore, to fully exploit the potential of this solu¬ 
tion, it requires the elaboration of a new hit position and TOE reconstruction 
methods 10 . Recently, one method of reconstruction in scintillator detectors 


based on the comparison of registered signals with respect to a library of syn¬ 
chronized model signals collected for a set of well-defined positions along the 
scintillator 0 is published. In this article we present similar method, however 
the comparison of measured signal is done with averaged model signals deter¬ 
mined as a function of position along the scintillator. This approach speeds up 
significantly the reconstruction in comparison to the previously used method. 


2 Library of synchronized model signals 


In order to create the library of model signals, we have used the setup com- 

scintilators with dimensions 300 x 19 x 5 mm^ and 
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posed of two BC-420 
Hamamatsu photomultipliers R49 98 |13| connected optically to their most dis¬ 
tant ends via optical gel EJ-550 |14| . General scheme of experimental setup 
used to build the library of signals is shown in Fig. 



Figure 1: Scheme of experimental setup where abbrevations PM and SDA denote 
Photomultilplier and Serial Data analyser (SDA6000A), respectively. 


The two detectors were irradiated with a collimated beam of annihilation 
quanta (FWHM ~ 1.5 mm) in steps of 3 mm. The ^^Na source was placed 
inside a lead collimator which could be moved along the scintillator by using a 
dedicated mechanical system. For each position a high statistics of signals with 
interval of 100 ps was collected by means of SDA. A coincident registration of 
signals from both detectors allowed for the photomultipliers noise suppression 
and the selection of annihilation gamma quanta. Exemplary sampled signal is 
shown in Fig. |^. To create a background-free library, first we have corrected 
all the signals for pedestal. For every signal the average value of voltage was 

^ Presently best TOP resolution was achieved with LSO crystals and amounts to about 
400 ps [^. 
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calculated in the noise region showed in Fig. by encircled red area. This 
computed average value was then used for pedestal correction for that particular 
signal. An exemplary signal after correction is presented in Fig. Dd. After 



(a) 



(b) 


Figure 2: (a) Measured signal from the database (b) The same signal after 
pedestal correction. 


pedestal correction, signals building the library were selected based on number 
of registered photoelectrons such that it is constructed from events with the high 
energy deposition. In Fig. we show an example of such spectrum obtained 
at central position. For every position we have cut all the signals for which 
the number of registered photoelectrons is lower than half of the number of 
photoelectrons corresponding to the Compton edge for 511 keV gamma quanta. 
Finally, database signals were synchronized by shifting their time scales in such 



Figure 3: Distribution of number of photoelectrons obtained from signals mea¬ 
sured irradiating the scintillator strip at central position(sohd black histogram). 
Dashed red histogram represents theoretical distribution fitted to the ex¬ 
perimental data using normalization and energy calibration constants as free 
parameters. Details of the fitting procedure are described in reference 
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a way that time of the gamma quantum hit inside the detector is the same for 
all events in the library. 

The shape of a model signal for each position from the database is determined by 
averaging pedestal corrected and synchronized signals. More details on model 
signals determination are given in the next section of this article. 
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3 Determination of the shape of model signals 

Function describing the approximate shape of the model signal was determined 
by averaging measured signals for a given hit position. Calculated average for 
signals registered simultaneously at both ends of the scintillator was treated as 
a reference in order to align all the database signals measured for a particular 
position, as it is shown in Figs. & and b. Such alignment is necessary to account 




Figure 4: (a) Example of a database signal before alignment to the average. 
Black curve represents the measured signal, while the average signal is shown 
by red dashed curve, (b) The same signal after rescaling. 


for the spread of signals in amplitude and time. To perform signal’s alignment in 
the database, we have defined for each measured signal at a given position a 
value. It was computed by comparing leading edge of the database and average 
signals. Comparison was performed by taking into account signal registered at 
the left and right side of scintillator simultaneously. The was calculated as 
a function of three parameters aL and aR\ 


X^{St,aL,aR) 


E jiModelLeftiVi) " tdbLeftijctLVi) " St))^ 
n. 


+E 


{tModelRightjVi) " ^db Right {{o^RVi) " 
m 


( 1 ) 


(5t is a time shift for sample points along time axis and o:r are normalisation 
factors for both signals (left and right) registered at both ends of scintillator. 
tModeiLeftiVi) and tModeiRight{Vi) denote time for model signal computed at 
left and right side for voltage Vi at their leading edge. idbLeftii^LVi) and 
tdbRight{{c^R^) Is the time computed for rescaled left and right signals at their 
leading edge, respectively. By minimization of the y^ value we obtained the best 
alignment between two compared signals and each database signal was rescaled 
using fit parameters o:r and St giving the lowest y^. Next, average of these 
rescaled signals was computed again leading to the model signal determination 
and the whole procedure was repeated until the changes were negligible. Ex¬ 
ample of model signals obtained at three different positions is shown in Eig. 
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Figure 5: Sample of model signals at three different positions: solid black line 
represents the left side of the scintillator, dashed red line center and dotted blue 
line the right side of scintillator. 

4 Validation and optimization of the time and 
position reconstruction 

In order to reconstruct hit position of gamma quanta in the detector the mea¬ 
sured signals are compared with each model signal from the library. This is done 
by minimizing value calculated at the leading edges of measured and model 
signals. Validation and optimization of this reconstruction method has been 
performed by utilizing the signals gathered at known positions with the same 
experimental setup which was used to build the library. In order to simulate the 
response of front-end electronics we have determined the times corresponding 
to predefined constant level thresholds [^: -60 mV, -120 mV, -180 mV and 
-240 mV. Analogously, we have chosen also four constant fraction thresholds 
to sample signals at: 1/8A^, 2/8A^, 3/84/8A^, where A^ is the amplitude 
of signal which is smallest among all four compared signals. To optimize the 
reconstruction we have considered two posibilities: 

1. when is a function of the time shift St only: 

4 

= y^X^ModelLeftiVi) — idbLeftiVi “ ( 2 ) 

4 

T ModelRight{y^i) 

i=l 

2. when x^ depends on St and normalisation factors ai, and aj^: 

4 

X^(St, aL, OiR) = y^X'^ModelLeftjVi) “ ^dbLeft{{oLLVi) - St)^ (3) 

i=l 

4 

+ '^^{tModelRightiVi) “ ^dbRight{{oLRVi) — St))‘^ . 
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tModeiLeft{Vi), t Model Right {Vi) and tdbLeft{Vi), tdbRight{Vi) denote times deter¬ 
mined at left and right side of the scintillator at threshold Vi for model and 
registered signals, respectively. The reconstructed hit position is the position of 
most similar signal from the library with respect to measured signal(i.e. model 
signal for which is minimaQ. 

The time of particle interaction is determined as a relative time between the 
measured signal and the most similar one from the library. This provides also 
determination of the gamma quantum time of flight {TOF) pT] : 

t firstStrip ^^firstStrip 

tsecondStrip ~ secondStrip (4) 

TOF ^secondStrip ^firstStrip i 

where 5t denotes shift in time for which the computed defined in Eq. and 
is lowest. An example of the y^ distribution calculated according to Eq. 



Eigure 6: Example of the y^ distribution defined in Eq. 

for one of signals sampled with constant levels and measured at central position 
of the strip is shown in Eig. One can see a clear minimum corresponding 
to z ^ 150 mm. Present version of reconstruction procedure do not take into 
account the measurement uncertainties. Therefore, the values are treated as 
arbitrary and errors of fitted parameters are determined directly from distribu¬ 
tions of differences between reconstructed and true values of time or position. 
Eigure. 1^ and b show distributions of differences between the true and recon¬ 
structed position using y^ defined by Eq. for constant fraction and constant 
levels discrimination, respectively. Resolution of position reconstruction is de¬ 
termined by fitting a gauss function to the presented distribution and obtained 
results give: cTz ~ 9.7± 0.3 mm for constant fraction sampling and cTz ~ 9.6 ± 
0.3 mm for constant levels. 


^ Since the degree of similarity is represented by the value. 
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Figure 7: Distribution of differences between the true and reconstructed position 
Az for signals measured at z = 150 mm for (a) constant fractions (b) constant 
levels. 
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Figure 8: Distribution of differences between the true and reconstructed position 
for signals measured at z = 150 mm for (a) constant fractions (b) constant 
levels 

In case when ^ is a function of (5t, and olr (Eq. obtained position 
resolution amounts to ~ 12.1 ± 0.4 mm and 10.9 ± 0.4 mm for constant frac¬ 
tion and constant level method, respectively. Corresponding lS.z distributions 
are shown in Fig. These resolutions were determined for signals measured 
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Figure 9: Position resolution as a function of the position of gamma quantum 
interaction for the case when (a) is a function of St only (b) is a function 
of Jt, aL and a^. 


at several positions along the scintillator as it is shown in Fig. The results 
indicate that this resolution does not change much with position. Results for 
TOF reconstruction for signals measured at z = 150 mm are shown in Fig. 
Corresponding resolutions are in this case equal to (Jtof ^ 103 ps for con¬ 
stant fractions sampling and (Jtof ^ 143 ps in case of constant levels. As it is 
shown in Fig. 11 results obtained for olr) amount to (Ttof ^ 132 ps 

(constant fractions) and (Trof ^ 119 ps (constant levels). 
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Figure 10: Distribution of differences between the true and reconstructed TOF 
for events registered at z = 150 mm. The reconstruction was done using ^ as a 
function of H with signal sampling at (a) constant fraction and (b) at constant 
levels 




Figure 11: Distribution of differences between the true and reconstructed TOF 
for signals measured at z = 150 mm. The reconstruction was done using as a 
function of 5i^aL^ and aR^ and with signal sampling at (a) constant fraction(b) 
constant levels. 
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In principle true value of TOF should be equal to zero when source was 
positioned in the middle between detection modules. However, due to different 
time offsets produced by electronics the reconstructed mean values of TOF may 
be different from zero. Again, we have performed studies of the TOF resolution 
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Figure 12: TOF resolution as a function of position along the scintillator: (a) 
= X^(<5i) and (b) = X^{St,aL,aR). 

as a function of position for both x^ functions as it is shown in Fig[l^ One can 
see that the determined resolutions is constant within the error bars over the 
full length of 30 cm long scintillator strip. 


5 Summary 


The preliminary results obtained during validation of the reconstruction method 
introduced in this article show that it is possible to obtain a spatial resolution 
of about 1.2 cm (a) for the gamma quanta hit position, and TOF resolution of 
about 125 ps (a). As regards the position resolution along the detector obtained 
result is few time worse than achieved by the commercial TOF-PET scanners. 
However as regards the TOF determination obtained result is better by about 
a factor of two with respect to resolutions achieved by the commercial TOF- 
PET tomographs characterized by typical field of views of about 16 cm and 
TOE resolution of about 230 ps (a) 18 . A further improvement is expected 


in the future by including measurement uncertainties and possible correlations 
between the times measured at different thresholds. 
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